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We systematically study the evolution of the Friedmann- Robertson-Walker (FRW) universe with 
a cosmological constant A and a perfect fluid that has the equation of state p = wp, where p and p 
denote, respectively, the pressure and energy density of the fluid, and w is an arbitrary real constant. 
Depending on the specific values of w, A and the curvature k of the 3-dimensional spatial space 
of the universe, we classify all the solutions into various cases. In each case the main properties 
of the evolution are studied in detail, including the periods of deceleration and/or acceleration, 
and the existence of big bang, big crunch, and big rip singularities. In some particular cases, the 
solutions reduce to those considered in some standard textbooks, where by some typos may be 
corrected. The methods used is simply the conservation law of kinetic and potential energies in 
classical mechanics, and undergraduate students can easily follow the analysis and apply it to the 
studies of other cosmological models of the universe. 

PACS numbers: 98.80.-k; 98.80. Jk;04.20.-q; 04.20.Jb 



I. INTRODUCTION 

Recent observations of supernova (SN) la reveal the 
striking discovery that our universe has lately been in its 
accelerated expansion phase !l|. Cross checks from the 
cosmic microwave background radiation and large scale 
structure all confirm this [2[ . Such an expansion was pre- 
dicted neither by the standard model of particle physics 
nor by the standard model of cosmology, and the under- 
lying physics still remains a complete mystery Q . Since 
the precise nature and origin of the acceleration have 
profound implications, understanding them is one of the 
major challenges of modern cosmology. As the Dark En- 
ergy Task Force (DETF) stated "Most experts believe 
that nothing short of a revolution in our understanding 
of fundamental physics will be required to achieve a full 
understanding of the cosmic acceleration." 

Within the framework of general relativity (GR), to 
account for such an acceleration, it requires the introduc- 
tion of either a tiny positive cosmological constant or an 
exotic component of matter that has a very large negative 
pressure and interacts with other components of matter 
weakly. This invisible component is usually dubbed as 
dark energy. For a perfect fluid with the equation of 
state, w = p/p, this implies w < —1/3, where p and p 
denote, respectively, the pressure and energy density of 
the fluid. On the other hand, a tiny positive cosmological 
constant is well consistent with all observations carried 
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out so far [5|. However, when we consider its physical 
origin, we run into other severe problems: (a) Its theo- 
retical expectation values exceed observational limits by 
120 orders of magnitude 7]. (b) Its corresponding energy 
density is comparable with that of matter only recently. 
Otherwise, galaxies would have not been formed. Consid- 
ering the fact that the energy density of matter depends 
on time, one has to explain why only now the two are in 
the same order, (c) Once the cosmological constant dom- 
inates the evolution of the universe, it dominates forever. 
An eternally accelerating universe seems not consistent 
with string/M-Theory, because it is endowed with a cos- 
mological event horizon that prevents the construction 
of a conventional S-matrix describing particle interaction 
Q. Other problems with an asymptotical de Sitter uni- 
verse in the future were explored in Q. 

In view of all the above, dramatically different models 
have been proposed, including quintessence [l(|, DGP 
branes [HI], and the f(R) models [Hj]. For details, see Q 
and references therein. However, it is fair to say that so 
far no convincing model has been constructed. 

To introduce such a fascinating subject to students, 
it is always challenging both physically and mathemat- 
ically. In this paper, our purposes are two-fold: (i) By 
studying the concrete models of the FRW universe, we 
shall show that such a subject can be easily understood 
in terms of classical mechanics. The method is very ba- 
sic, and undergraduate students with some knowledge 
of classical mechanics can easily follow. As a matter of 
fact, the only knowledge required is the conservation law 
of total energy of a classical particle with mass m moving 
under a potential V(x) [Hj], 

^mi 2 + V(x) = E, (1.1) 

where x = dx/dt and E is the total energy of the system, 
which is conserved without external force. Then, taking 



2 



derivative of the above equation with respect to t, we find 
that 



dV(x) 
dx 



(1.2) 



Thus, once the potential V(x) is known in terms of x, 
one can immediately tell if the particle is accelerating or 
decelerating, without integrating Eq. (|l.l|) explicitly. In 
addition, one can also easily determine the range of x 
that the motion allows. Therefore, if the problem of the 
evolution of the universe can be expressed in the above 
form, we can use the methods of classical mechanics to 
study its evolution, and classify all the possible models 
of the universe, (ii) The second purpose of this paper is 
to correct some typos appearing in some standard text- 
books. 

The rest of the paper is organized as follows: In Sec. 
II, we consider the Friedmann equation coupled with a 
cosmological constant and a perfect fluid with the equa- 
tion of state p = wp for any given curvature k. After 
writing it in the form of Eq. , we study the potential 
V{x) case by case, and deduce the main properties of 
each model of the universe. In Sec. Ill, we present our 
main conclusions. 

It should be noted that classification of a (non- 
relativistic) matter coupled with a dark energy was con- 
sidered recently in [3], and the corresponding Penrose 
diagrams were also presented. 

In this paper, we shall use the notations and conven- 
tions defined in [l5[. 



II. CLASSIFICATION OF THE FRW UNIVERSE 

The FRW universe is described by the metric [HI], 
dr 2 



ds 2 



dt 2 -a 2 {t) 
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in the spherically symmetric coordinates x a = 
{t,r,6,(j)} 7 (a = 0,1,2,3), where k denotes the curva- 
ture of the three-dimensional spatial space of constant 
t, and can be set k = 0,±1, without loss of generality. 
a(t) is the expansion factor of the universe. It should be 
noted that Eq. (|2.1[) is invariant under the translation, 



t' = t - t. 



(2.2) 



where t s is a constant. In the following we shall use this 
gauge freedom to fix the origin of the timelike coordinate 
t. The expansion factor a(t) of the universe is determined 
through the Einstein field equations 



1 2 

Rab — -^R9ab = K T a b + ^9ab 



(2.3) 



where k 2 [= 8ttG/c 4 ] is the Einstein coupling constant, A 
denotes the cosmological constant, and T a b the energy- 
momentum tensor of the matter field(s) filled in the uni- 
verse. For a perfect fluid, we have 



Tab = (P+P) UaU b ~ P9ab, 



(2.4) 



where u a = de notes the four-velocity of the fluid. It 
can be shown [la] that the Einstein field equations (|2.3p 
for the metric (|2.1|) and energy-momentum tensor (|2.4|) . 
have only two independent components, which can be 
cast in the form, 



H 2 = 



a 
a 



8irG 
4nG 
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f 3 A ^' 

(P + 3 P ) 



1 



-A, 



(2.5) 
(2.6) 



where H = a/a. Note that in writing the above equation, 
we had chosen units such that the speed of light is one. 
On the other hand, the conservation law of matter fields, 
V a T a6 = 0, yields 



p + 3H (p + p) = 0. 



(2.7) 



It can be shown that this equation is not independent, 
and can be obtained from Eqs. (|2.5[ ) and (|2.6|) . 

Note that we have three unknowns, a, p and p, but 
only two independent equations. Thus, to close the sys- 
tem, we need to have one more equation. Usually, this is 
given by the equation of state of the matter field. In this 
paper, we shall consider the case where 



(2.8) 



p = wp, 



where w is an arbitrary real constant. Inserting Eq. (|2.8[) 
into Eq. (|2.7p , we find that it allows to integrate once and 
gives, 



(2.9) 



where po and ao are the integration constants. Since 
Po represents the energy density when a = ao, we shall 
assume that it is strictly positive po > 0. Without loss 
of generality, we can always set ao = 1. Then, it can be 
shown that the Friedmann equation (|2.5p can be cast in 
the form of Eq. (fi~T|) with m = 1, E = 0, x(t) = a(t), 



where 



with 



V(a) 



-d 2 + V(a) = 0, 



-k - -Aa 2 
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(2.10) 



(2.11) 



(2.12) 



When w — the problem reduces to the one treated in 
[l5l ]. To study the problem further, we consider the cases 
k = 0, ±1 separately. 
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A. k = 
When k = 0, Eg. (|2JTj) reduces to 



7(a) 



-A^ 



C 



,l+3w ' 



(2.13) 



It is found convenient to consider the cases where A > 
0, A = and A < separately. 



1. A > 

When A > 0, Eq. (|2.13p can be written as 



V(a) 



Aa 2 1 



C 



,1+3™ ) ' 



(2.14) 



where C = 6C/|A| > 0. It is found that, depending on 
the value of w, the evolution of the university can be 
significantly different. Thus, we shall further distinguish 
the following sub-cases: 

(i)w>-~, (ii)w = ~, (iii)-Kw<~, 



(iv) w = — 1, (v) w < — 1. 



(2.15) 



Case A. 1.1) w > — |: Then, we find that V(a) 
is strictly negative, and V{a) — > — oo for both a = 
and a — > oo. It also has a maximum at a = a m = 
(3 (1 + 3w)C/A) 1/m+w)) , for which 



< 0, a < a.„ 
= 0, a = o„ 
> 0, a > a v 



(2.16) 



Fig. [T] schematically shows the potential. Therefore, 
in this case, the evolution of the universe is dominated 
by matter in the early time and a(t) scales like a(t) oc 
/;2/[3(i+tu)]^ f or ^^jjiojj a < 0. As the universe expands to 
a = <hni it researches the turning point, after which it 
starts to expand acceleratingly, i.e., d(a > a m ) > 0. The 

universe is asymptotically de Sitter, a(t) oc eV A / 3 *. 



Case A. 1.2) w 



V(a) 



In this case, we find 



--Aa 2 -C 
6 

— C, a - 
—oo, a - 

dV(a) 
da 



oo, 



-Aa > 0. 



(2.17) 



The middle line of Fig. [T] is the potential for this case. 
Thus, in this case the universe is always accelerating ex- 
cept for the initial moment a = 0. On the other hand, 



from Eq. 



we find that p oc a 



oo as a — > 0, that 



is, a big bang singularity still occurs at a = 0. 




FIG. 1: The potential given by Eq. (f2TT4|) for k = and A > 0. 



Case A. 1.3) -1 < w < -| 

the potential is that given by Fig. [TJ and 



In this case, we find that 



dV(a) 
da 



Aa + (3\w\ - 1) Ca 3 H- 2 > o, (2.18) 



where the equality holds only at a — for w = —2/3. 
Thus, in this case the universe is always accelerating. 
Note that in the present case a = still represents a big 
bang singularity, as can be seen from Eq. (|2.9[) . It is also 
interesting to note that at a — we have a = d = 0, 
except for w = —2/3. Thus, when w ^ —2/3 the point 
a = is a stationary point. However, it is not stable, 
and any perturbations will make the universe to expand. 
When w = — 2/3 we have d(a — 0) = and d(a = 0) = 
(3|k;| — 1) C > 0. Therefore, in the latter case a = is not 
stationary, and the positive force will lead the universe 
automatically to expand. 

Case A. 1.4) w = — 1: In this case, the matter acts as 
a vacuum energy, and the potential is given by 



V(a) 



A eff a 2 , 



(2.19) 



where A e ff = A + 6C. Therefore, in this case the universe 
is de Sitter, and 



o(t) 



V A "/// 3 (*- to ) 



(2.20) 



Recall that the de Sitter space is free of any kind space- 
time singularities at a — as well as at a = oo [To] ]. 

Case A. 1.5) w < —1: In this case, the behavior of 
the potential V(a) and a(t) are similar to the case — 1 < 
w < —1/3, except for the fact that now p oc a 3 *^' 1 "' -1 '' is 
not singular at a = 0, although it is at a — oo, which 
is usually called a big rip singularity. At a = we have 



FIG. 2: The expansion factor a(t), the acceleration d(t), and 
the energy density p(t) for k = and A > 0. The spacetime 
has a big bang singularity at t = for w > — 1. It is de 
Sitter for w = — 1, which is free of any kind of spacetime 
singularities. When w < — 1, a big rip singularity occurs at 
t = t 3 , at which we have a (t s ) — oo = p (t s ). 



a = a — 0, that is, in this case this point also represents 
a unstable stationary point. 

Finally, we note that for all the cases with k = the 
corresponding Friedmann equation can be integrated ex- 
plicitly, and the correspond solutions are given by, 



'6C\ 1/2 
a(t) = {( — ) sink 



x /3A , 
(l+w)J—(t-t s ) 



for where t s is given by 

i rr if ax n 



(2.21) 



(2.22) 



so that a(t = to) = 1. For w > — 1, without loss of gen- 
erality, we can always use the gauge freedom of Eq. (|2.2p 
to set t s = 0, so that the big bang singularity occurs 
at t = 0. This will be the case in the rest of this pa- 
per. When w — —1, the solution is that of de Sitter, 
given by Eq. (12.20|) , which is free of any kind of spacetime 
singularities, and the solution is valid for any t, that is, 
t € (—00,00). For the cases where w < —1, we shall 
not do such a translation, so that the big rip singularity 
happens exactly at t = t s . Note that when w < — 1, from 
Eq. (|2.21[ ) we find that the solution is valid only when 
t G (— 00, ts), for which we have 



»(*) = 




1) {ts - t) 



3(1" 



(w < -1) 



(2.23) 



Since now we have p oc a 3 ^' 1 "' -1 ^ we find that the space- 
time is indeed not singular at a(t = —00) = 0, but sin- 
gular at a(t — t s ) = 00. 

In Fig. [2] we summarize the main properties of the 
solutions for k — and A > with different w. 




FIG. 3: The potential given by Eq. (|2~T4|) for k = and A = 0. 



2. A = 



In this case, we have 
V(a) = 



C 



,l+3iu 



<o, 



-(3tu+l)- 



C 



j2+3w 

W > -1/3, 
w = -1/3, 
w < -1/3. 



(2.24) 



Fig. [3] shows the potential V(a). Similar to the last case, 
now we can also integrate the Friedmann equation (|2.5p 
to obtain the explicit solutions of a(t) and p(t), 



a(t) = 

P(t) = 
where po = 



3(l + w)Jf(t-i,) 



/2C(t-t ) 



W 7^ 
10 = 



Po 



w 7^ 



2p /[9{l + w) 2 C], and 



in 



1 



(2.25) 



(2.26) 



When w < —1, Eq. (|2.25p shows that, to keep a(t) real 
and positive, we must require t G (— oo,t s ). As £ — ► 
—00, both a(t) and p(t) vanish, while when t — > t s all of 
them become unbounded, that is, a big rip singularity is 
developed there. 

In Fig. 2] we summarize the main properties of the 
solutions for k — and A = with different values of w. 
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FIG. 4: The expansion factor a(t), the acceleration a(i), and 
the energy density p(t) for k — and A = 0. There is a 
big bang singularity at a = for all the cases with w > — 1. 
The spacetime is de Sitter for u? = — 1. When to < — 1, a 
big rip singularity is developed at t = t s , at which we have 
a(t s ) = p(i s ) = oo. 



3. A < 



In this case, we have 



V(a) 



\A\a 2 



,l+3tu : 



which has been shown in Fig. [5] for various values of w. 
It can also be shown that the corresponding expansion 
factor and the energy density are given, respectively, by 




a(t) = 
P(t) = 

for w ^ — 1, and 



-Y(l + w) (t-Q 



Po 



+ (t-t s ) 



(2.28) 



a(t) = e VW3(*-*o), 
p(t) - po, (2.29) 
for w = —1, where A e y/ = |A| — 6C > 0, but now with 



l + w 




Po 



\Mpq 

6C ' 



(2.30) 



As what we did in the previous cases, using the gauge 
freedom (|2.2|l we shall set t s = for w < — 1, while keep 
it as it is for u> < — 1. 

Case A. 3.1) w > — |: In this shown in Fig. 

we have V (a > a m ) > 0, so the motion of a > a m is 
forbidden. If the universe starts to expand from the big 
bang where a = 0, when it expands to its maximal radius 
a m it will starts to collapse until a big crunch singularity 
is developed at t = 2t m , as shown in Figs. [S] and [51 where 
dm = a (tm)i and is given by 



6C \ 3 
|A|. 



(2.31) 




(2-27) FIG. 5: The potential given by Eq.([22ZJ) for k = and A < 0. 



During the whole process, the universe is always deceler- 
ating, 



dV(a) 

a = < 0, 

da 

as can be seen from Fig. [SJ 

Case A. 3. 2) w = — |: In this case, we find 



(2.32) 



V(a) 



(i 



\A\a 2 -C 



—C, a = 0, 
oo, a — > oo, 



-_Li=-_ A o<0. 
da 3 



(2.33) 



Therefore, similar to the last case, the universe expands 
from the big bang singularity at a = until its maximal 
radius a m , given by Eq. (|2.31[) with w — —1/3, and then 
starts to collapse until a big crunch singularity is formed 

Case A. 3. 3) — 1 < w < — |: In this case, from Fig. O 
we can see that V(a = 0) = = V (a m ), and the motion 
is also restricted to a < a m . However, there is a fun- 
damental difference between this case and the last two 
cases: The potential V(a) has a minimum at a = a m i n , 
at which we have dV (a m , n ) /da — 0. The universe is 
initially accelerating. But, when it expands to a = a m in, 
it starts to decelerate until a = a m , at which its expand- 
ing velocity becomes zero, and afterwards it will starts 
to collapse, until a big crunch singularity is developed at 
a = 0, as shown in Fig. [S] 

Case A. 3. 4) w = — 1: In this case, we have 



(2.34) 
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FIG. 6: The expansion factor a(t), the acceleration d(t), and 
the energy density p(t) for fc = and A < 0. The spacetime 
is singular at a = for all the cases with w > — 1, (a big bang 
singularity), it is de Sitter for w = — 1. When iu < —1, a 
big rip singularity is developed at t = t s , at which we have 
a(i s ) = p(£ s ) = oo. 



Therefore, now there is solution only when A e ff > 0, for 
which the universe is dc Sitter, and 



a(t) 



= P V A =/// 3 ( t -*o) 



(2.35) 



where A eff = |A| - 6C. 

Case A. 3. 5) w < — 1: In this case, there is a minimum 
Q-rnm for which V (a < a m in) > 0. Therefore, in contrast 
to the previous case, now the motion of the universe is 
restricted to a > a m i n . The universe starts to expand 
from a — a m i n until a (t s ) = oo within a finite time, 
whereby a big rip singularity is formed, as shown in Fig. 

El 



B. k = 1 



In this case, the potential given by Eq. (|2.11j ) can be 
written as 



V(a) 



1 1 

2 ~ 6 J 



-Aa 2 




(2.36) 



3(l+u>) 



Following the case of k = 0, 



where C = (12C/A) 
we also distinguish the three cases, A > 0, A = and 
A < 0. 



1. A > 

When A > 0, it is found convenient to further divide 
it into the five sub-cases defined as in Eq. (|2.15| ). 

Case B.l.l) w > — |: In this case, it can be shown 
that for any given w and po there always exists a critical 
value A c and radius a m that satisfy the conditions, 

V(a m ,w,p , A c ) = 0, V (a m , w,p , A c ) = 0, (2.37) 

where a prime denotes the ordinary differentiation with 
respect to a. It can be shown that the solutions of the 




FIG. 7: The potential given by Eq. (p36]) for k = 1, 
and A > 0, where A c = A c (w, po). 



above conditions are 



w > 



A, 



= [3{l + w)C] 
'l + 3w 



1 + w 



[3{i + w )cy 



(2.38) 



As to be shown below, the solutions with A > A c have 
quite different properties from the ones with A < A c . 
Therefore, in the following we shall further distinguish 
the three different cases, A > A c , A = A c and A < A c . 

Case B. 1.1.1) w > — |, A > A c : In this case, the po- 
tential V(a) is always negative for any given a, as shown 
in Fig. [7] Therefore, the corresponding solutions have no 
turning point. If the universe initially starts to expand 
from a big bang singularity at a = 0, it will expand for- 
ever, as shown by Fig. [8] However, the potential has a 
maximum at a = a m , for which we have 



< 0, a < a m , 
= 0, a = a m , 
> 0, a > a m , 



(2.39) 



that is, the universe is initially decelerating. Once it 
expands to a = a m , it starts to expand acceleratingly. 

Case B.l.l. 2) w > — |, A = A c : In this case, there 
exists a static point a m , at which we have V (a m ) — 
V (a m ) — 0, as one can see from Fig. [7] Therefore, if 
the universe starts to expand from the big bang at a = 0, 
it will expand until a = a m . The universe is decelerating 
during this period. Since at the point a = a m , we have 
a = = a, the universe will become static, once it reaches 
this point. However, it is not stable, and with small 
perturbations, the universe will either collapse until a 
singularity is developed at a = or expand forever with 
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FIG. 8: The expansion factor a(t), the acceleration a{i), and 
the energy density p(t) for k = 1, w > — | and A > 0, where 
A c is given by Eq. (|2.38[) . A big bang singularity occurs at 
a — in all cases with A > A c . In the first sub-case of 
OA < A c , both big bang and big crunch singularities occur, 
while in the second sub-case the spacetime is free of any kind 
of spacetime singularities. 



a > 0. If the universe initially at a — ai > a m , from Fig. 
[7] we can see that it will expand forever. Since V'(a) is 
always negative, so the universe in this region is always 
accelerating. 

Case B.l.1.3) w > — |, < A < A c : In this case, 
V{a) = has two real and positive roots, a\ and a-i, 
as shown in Fig. [7J Without loss of generality, we as- 
sume that d2 > a%. Since V(a) > for a G (01,02), the 
motion is forbidden in this region. Similar to the last 
case, depending on the initial conditions, the universe 
can have quite different evolutions. In particular, if it 
starts to expand from the big bang singularity at a = 0, 
it will expand until a — a\, at which we have a — and 
d < 0. Since d < at this point, the universe will start 
to collapse afterwards, until a big crunch singularity is 
developed at a = 0. If the universe starts to expand at 
a-i > 0-2, it will expand forever. In the latter case, the 
universe is always accelerating, as can be seen from Fig. 



7. 



Case B.1.2) 



-|: In this case, we have 




FIG. 9: The potential given by Eq. lpOOj) for k = 1, w < -§ 
and A > 0. 
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FIG. 10: The expansion factor a(t), the acceleration d(t), 
and the energy density p(t) for k = 1, w < — | and A > 0. 
A big bang singularity occurs only in the case w = —1/3 and 
C > 1/2. In the case w < — 1, a big rip singularity occurs at 
a = 00. 



V(a) = -(1-C) 



G 



-Aa 2 



(2.40) 



Thus, depending on the value of C (po), the motion 
of the universe will be different. In particular, when 
C (po) < 1/2, there exists a minimal a m i n , for which 
V (a < a m in) > 0, that is, the motion in the region 
< a < a m i n is forbidden, as shown in Fig. [5J When 
C (po) ^ 1/2, the universe can start to expand from the 
big bang singularity at a = 0. In all the cases we have 
V'(a) < 0, so that the universe is always accelerating [cf. 
Fig. Q3J|. 

Case B.1.3) — 1 < w < — 4: In this case, we find that 
V'{a) is strictly negative for any a > with V(0) = 1/2. 
Therefore, similar to the case w = —1/3 and C < 1/2, 
there exists a minimal a m i n , for which V (a < a m i n ) > 0, 
and the motion in the region < a < a m i n is forbidden, 
as shown in Fig. [9l Thus, in the present case the universe 
starts to expand from a radius > a min until a = oo 
without turning point. Again, because now V'(a) < for 



any a > a m i n , the universe is always accelerating. How- 
ever, there is no any kind of singularities to be developed, 
either big bang, big crunch, or big rip, as shown by Fig. 

Case B.1.4) w = — 1: In this case, the potential is a 
simple parabola, 



V{a) 



(A + 6C) a 2 



(2.41) 



schematically shown by the top curve in Fig. [5] As a 
result, the motion is similar to the last case, except for 
the fact that now p = pa. 

Case B.1.5) w < —1: In this case, we also have 
V'{a) < and there exists a finite radius, a m i n , such that 
when a < a m i n we have V(a) > 0, and when a > a m i n we 
have V(a) < 0. The only difference is that in the present 
case there is a big rip singularity that happens at a = oo, 
as now we have shown in Fig. [TU1 
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FIG. 11: The potential given by Eq. lfPS]) for fc = 1 and 
A = 0. 



A = 

In this case, we find that 



V (a) = ^ - (2-42) 

Fig. [TT] shows the potential for various values of w, from 
which we can see that when w > —1/3, the motion of the 
universe is restricted to a < a m , where a m is the solution 
of V(a) — 0. The universe starts to expand at the big 
bang singularity a = until the turning point a = a m . 
Afterwards, it will start to collapse until a big crunch 
singularity is developed at a — 0, as shown in Fig. [T2l 

When w = —1/3, there is motion only for C > 1/2, 
for which the universe expands linearly from a big bang 
singularity at a — with a = 0. 

When —1 < w < —1/3, the motion is possible only 
for a > a m i n , as shown in Figs. llll and [T2"l The universe 
starts to expand from the initial point a, > a m i n with a > 
0. No turning point exists, so the universe will expand 
forever. During the whole process, the matter density 
remains finite, so no singularity exists in this case. 

When w < —1, it can be shown that the motion for 
a < CLmin is also forbidden. As a result, no big bang 
singularity exists in the present case. But, a big rip sin- 
gularity will be developed as a — > oo, as shown by Fig. 
IT2"1 In the whole process, we have a > 0. 
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FIG. 12: The expansion factor a(t), the acceleration a(t), and 
the energy density p(t) for k = 1 and A = 0. There are both 
big bang and big crunch singularities in the case w > —1/3, 
while only a big bang singularity occurs in the case w — —1/3. 
There is no singularity in the cases with — 1 < w < —1/3. A 
big rip singularity occurs at a — oo for w < — 1. 

3. A < 

In this case, we have 

V{a) = \ + \\k\a*--^, (2.43) 

which has the properties as shown by Fig. [T31 In partic- 
ular, when w > —1/3, we find that the universe starts to 
expand from a big bang singularity at a = until a max- 
imal radius a m where V (a m ) = 0. Afterwards, the uni- 
verse starts to collapse, and finally a big crunch is devel- 
oped at t = 2t m where t m is determined by a rn = a (t m ). 
In the whole process, the universe is decelerating, as 
shown by Fig. RBI 

When w = —1/3, the potential is non-negative for C < 
1/2, so the motion is forbidden. When C > 1/2 we have 
V(a) < for a < a m , where a m is the root of V(a) = 0, 
as shown in Fig. [T31 Thus, the motion now is possible 
in the region a < a rn , for which the universe starts to 
expand from a big bang singularity at a = 0. Once it 
reaches its maximal at a m , it starts to collapse until a 
big crunch is developed at t = 2t m . 

When —1 < w < —1/3, similar to the last case, the 
potential is negative only for a < a m , as shown in Fig. 
1131 In particular, a big bang (crunch) singularity happens 
at t = (t = 2t m ). The difference is that now there exists 
a time t max so that for < t < t max or 2t m - t max <t< 
universe is accelerating, while during the time 
t m -t ma x <t < 2t m -t max it is decelerating, where t max 
is the root V (t max ) = 0. 

When w = — 1, the potential is non-positive only for 
C > |A|/6 and a > a m i n , where a m i n is the root of 
V(a) = 0, as shown in Fig. [131 Therefore, in this case 
the universe starts to expand always at an initial radius 
Ori > a m in- The universe will expand forever with a > 0. 
However, the spacetime is not singular even when a = oo. 

When w < — 1, the potential is non-positive only for 
a > a-min, where again a m i n is the root of V(a) = 0, as 
shown in Fig. 1131 The evolution of the universe in this 
case is similar to the last one, except for that now a big 
rip singularity will be developed at a — oo. 



A V(a) 





FIG. 15: The potential given by Eq. (|27H)l for k = -1 and 
A > 0. 



FIG. 13: The potential given by Eq. (p^3|) for k = 1 and 

A = 0: (a) for to > -1/3; (b) for'™ = -1/3 and C > 1/2; (c) 
for -1 < w < -1/3; (d) for w = -1 and C > |A|/6; and (e) 
for ui < — 1. 
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In this case, we find that V(a) — > — oo as a — > oo, and 
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FIG. 14: The expansion factor a{t), the acceleration a{t), and 
the energy density p(t) for k = 1 and A < 0. There are both 
big bang and big crunch singularities in the case with w > — 1, 
while only a big bang singularity occurs in the case w = — 1, 
and no singularities for the case w = — 1, while there is a big 
rip singularity at a = oo for w < — 1. 



C. k = -l 



f -1/2, w<-l/3, 
V(o) = ^ -(1/2 + C), io = -1/3, (2.45) 
t — oo, w > —1/3, 

when a — > 0, as shown by Fig. [15] Thus, when w < —1/3, 
the potential has a maximum at a m where V (a m ) = 0. 
The universe starts to expand from a big bang at a = 0. 
Initially, it is decelerating, a < 0. However, when it 
expands to a m , it turns to expand acceleratingly, a > 0, 
as shown in Fig. [TBI When — 1 < w < —1/3, the universe 
expands from a big bang at a — until a — oo, and 
there is no turning point. It expands with a > in the 
whole process. The case of w = — 1 is similar to the case 
of —1 < w < —1/3, except that the spacetime is not 
singular either at a = or at a = oo, as shown in Fig llGI 
When w < — 1, one can see that the universe starts to 
expand from a = with a > for any given a. There 
is no big bang singular at a = 0, but there is a big rip 
singularity at a = oo. 



When k = — 1, the potential is given by 

V(a) = — - Ika 2 - (2.44) 

To study the motion of the universe in this case, it is 
also found convenient to distinguish the three cases A > 
0, A = and A < 0, and in each case there are five 
sub-cases with different choices of w. 



2. A = 

In this case, we have 

V(a) = -I- -fx;- < 0, 

a = -^+1£= 0, . = -1/3,(2.46) 
a l>0, w<-l/3, 
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FIG. 16: The expansion factor a{t), the acceleration a{t), and 
the energy density p(t) for fc = — 1 and A > 0. There are a 
big bang singularity for w > —1, no singularity for w = — 1, 
and a big rip singularity at a — oo for w < — 1. 




-1/2-C 



FIG. 17: The potential given by Eq. flggg} for fc = -1 and 
A = 0. 



when a € [0, oo), as shown by Fig. [T7J We also have 



p(a) 



as a — * 0, and 



oo, 


tu > 


-1/3, 




Po, 


w = 


-1/3, 


(2.47) 


o, 


w < 


-1/3, 




o, 


w > 


-1/3, 




Po, 


w 


-1/3, 


(2.48) 


oo, 


V < 


-1/3, 





as a — > oo. Fig. 1181 shows the motion of the universe for 
each given w. 



3. A < 



In this case, we have 
Via)- 



1 1 

2 + 6 



|A|a 2 



C 



,1+3mi ' 



(2.49) 
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FIG. 18: The expansion factor a(t), the acceleration a(t), and 
the energy density p(t) for k = —1 and A = 0. There are a 
big bang singularity for w > — 1, no singularity for w = — 1, 
and a big rip singularity at a — oo for w < — 1. 



Depending on the vales A, C and w, the potential will 
have quite different properties. In the following we shall 
study all of them case by case. 

Case C.3.1) w > —1/3: In this case, the potential 
is shown schematically in Fig. I19[ from which we can 
see that it is non-positive only for a < a m , where a m 
is the positive root of V(a) — 0. Clearly, in this case 
there is a big bang singularity at a = 0, from which the 
universe starts to expand until a = a m . Afterwards, it 
will collapse so that finally a big crunch singularity is 
developed at t = 2i m , at which we have a (2t m ) = 0, as 
shown by Fig. 

Case C.3.2) w = —1/3: The potential in this case 
is similar to the last case, except for the fact that now 
1/(0) = -1/2 - C, as shown in Fig. [H The motion of 
the universe is qualitatively the same as that in the last 
case, as shown by Fig. [501 

Case C.3.3) — 1 < w < —1/3: In this case the po- 
tential has a minimum at a = a m i ni as shown in Fig. 1191 
for which we find that a < for a < a m i n , and 'a > for 
a > Cbmini a s shown by Fig. [2D] 

Case C.3.4) w = —1: In this case, depending on 
the ratio 6C/|A|, there are three distinguished sub-cases. 
When 6C/|A| < 1, the potential is non-positive only when 
a < a m where a m = [3/(|A| — 6C)] 1 / 2 , as shown by Fig. 
1211 Then, the universe starts to expand from a — until 
a = a m . Afterwards, it will start to collapse until a = 
again. But in the whole process, no spacetime singularity 
is developed. So, the universe is oscillating between a = 
and shown in Fig. 1221 

When 6C/|A| = 1, we find that V(a) = -1/2, and the 
universe expands linearly starting from a — 0. There is 
no turning point, and no spacetime singularity, as shown 
by Figs. EH and H 

When 6C/|A| > 1, we find that V(a) < -1/2 for any 
given a. Then, starting from a = 0, the universe expands 
always acceleratingly (a > 0) until a = oo, as shown by 
Figs. |2"T1 and |2"2"1 No spacetime singularity is developed 
during the whole process. 

Case C.3.5) w < —1: In this case, it can be shown 
that for any given po there always exists a critical value 
A c so that V{a) = has two root positive roots when 
|A| > |A C |; one positive root when |A| = |A C |; and no 
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FIG. 19: The potential given by Eq. (pT49"j) for k = — 1, A < 
and w > —1. 
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FIG. 20: The expansion factor a(t), the acceleration a{t), and 
the energy density p(t) for fc = — 1, A < and w > — 1. There 
are both big bang and big crunch singularities for all the cases 
with w > — 1. 

positive root when |A| < |A C |, as shown by FigfT^l where 
A c is the solution of the equations V(a m ,A c ) — 0, and 
V (a m , A c ) = 0. It can be shown that it is given by, 

Ac = (Hrf ) I 3 (M - 1) . (2.50) 

Case C.3.5.1) |A| > |A C |: In this case, the potential 
is positive in the region a\ < a < a 2 , where ai and a 2 
are the two positive roots of V(a) — with a 2 > a\. 
Therefore, the motion of the universe now is restricted 
to the regions < a < a\ and a > a 2 , depending on 
its initial condition. If the universe starts to expand at 
a = 0, it will expand until its maximal radius a = a±, and 
then collapse until a = 0. In the whole process, we have 
a < 0. Since at a = the spacetime is not singular, so 
the universe will start to expand again. This process will 
be repeating endlessly, as shown in Fig. [M] However, 
if it starts to expand at a radius dj > a 2 , the universe 




FIG. 21: The potential given by Eq. (pT49j) for k = — 1, A < 
and w — — 1. 
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FIG. 22: The expansion factor a(t), the acceleration a(t), and 
the energy density p(t) for k — — 1, A < and li) = — 1. The 
spacetime is not singular in any of these cases. 

will expand forever and never stops, as now a > for 
any given a > a 2 . A big rip singularity will be finally 
developed at a = oo, since now we have p — > oo, as 
a —>■ oo. 

Case C.3.5.2) |A| = |A C |: In this case, there ex- 
ists a static point a m , at which we have V (a m ) — 
V (a m ) = 0, as one can see from Fig. [23l where 
a m = [3{\w\ -l)C] _1/(3|w|_1) . Therefore, if the uni- 
verse starts to expand from a = 0, it will expand until 
a = a m with a < 0. Since at the point a = a m , we have 
a = = a, the universe will become static at this point. 
However, it is not stable, and with small perturbations, 
the universe will either collapse until a = or expand 
forever with a > 0. It should be noted that the space- 
time is not singular at a = 0. So, if it collapses, it will 
start to expand again when the point a = is reached. 
If the universe initially at a — ai > a m , from Fig. [55] we 
can see that it will expand forever. Since V [a) is always 
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FIG. 23: The potential given by Eq. (p^9|) for k = -1, w < 
— 1 and A < 0, where A c is given by Eq. (|2.50[) . where (a) is 
for |A| > |A C |; (b) is for |A| = |A C |; and (c) is for |A| < |A C |. 
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FIG. 24: The expansion factor a{i), the acceleration a(t), and 
the energy density p(t) for k = — 1, w < — 1 and A < 0, where 
A c is given by Eq. (|2.50[) . There are big rip singularities in all 
the cases, except for the first sub-cases of |A| > |A C | and 
|A| = |A C |. 

negative, so the universe in this case is always accelerat- 
ing. A big rip singularity will be also finally developed 



at a — oo. 

Case C.3.5.3) |A| < |A C |: In this case, the potential 
V(a) is always negative for any given a, as shown in Fig. 
1231 Therefore, the corresponding solutions have no turn- 
ing point. If the universe initially starts to expand from 
a = 0, it will expand forever. However, the potential has 
a maximum at a — a m , for which we have 

r < 0, a < a m , 
a = < = 0, a = a m , (2-51) 
{ > 0, a > a m , 

that is, the universe is initially decelerating. Once it 
expands to a = a m , it starts to expand acceleratingly. 
Similar to the last two cases, a big rip singularity will be 
also finally developed at a = oo. 

III. CONCLUSIONS 

In this paper, we have systematically studied the so- 
lutions of the Friedmann-Robertson- Walker (FRW) uni- 
verse with a cosmological constant and a perfect fluid 
that has the equation of state p — wp, where p and p 
denote, respectively, the pressure and energy density of 
the fluid, and w is an arbitrary real constant. Writing 
the motion of the universe in the form of Eq. (|l.l[) ,we 
have been able to classify all the solutions according to 
the different values of k, po, w and A, by simply us- 
ing the knowledge of one-dimensional motion in classical 
mechanics [li|, where k[= 0, ±1] denotes the curvature 
of the FRW 3-dimensional spatial space, po the energy 
density of the matter field when a = 1, and A the cos- 
mological constant. All these solutions are classified and 
presented in Figs. 0HllSiai^II21IiaiI3IIBlEDIE2l 
and 1241 Some particular cases were already discussed in 
various standard textbooks, whereby some typos may be 
corrected. 

In this paper, we have studied the FRW universe in 
the Einstein theory of gravity. Clearly, such studies can 
be easily applied to all the models where the evolution of 
the universe can be cast in the form of Eq. ljl.ip [3], in- 
cluding brane worlds in string/M theory [16j |. asymmet- 
ric branes [l3], and cosmological models in the Horava- 
Lifshitz gravity [HI]. 
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